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ABSTRACT 

We have undertaken the study of the elemental abundances and kinematic properties of a metal- 
poor sample of candidate thick-disk stars selected from the RAVE spectroscopic survey of bright stars 
to differentiate among the present scenarios of the formation of the thick disk. In this paper, we report 
on a sample of 214 red giant branch, 31 red clump/horizontal branch, and 74 main-sequence/sub-giant 
branch metal-poor stars, which serves to augment our previous sample of only giant stars. We find 
that the thick disk [a/Fe] ratios are enhanced, and have little variation (< 0.1 dex), in agreement 
with our previous study. The augmented sample further allows, for the first time, investigation of the 
gradients in the metal-poor thick disk. For stars with [Fe/H] < —1.2, the thick disk shows very small 
gradients, < 0.03 ± 0.02 dex kpc -1 , in a-enhancement, while we find a +0.01 ± 0.04 dex kpc -1 radial 
gradient and a —0.09 ± 0.05 dex kpc -1 vertical gradient in iron abundance. In addition, we show that 
the peak of the distribution of orbital eccentricities for our sample agrees better with models in which 
the stars that comprise the thick disk were formed primarily in the Galaxy, with direct accretion of 
stars contributing little. Our results thus disfavor direct accretion of stars from dwarf galaxies into 
the thick disk as a major contributor to the thick disk population, but cannot discriminate between 
alternative models for the thick disk, such as those that invoke high-redshift (gas-rich) mergers, heating 
of a pre-existing thin stellar disk by a minor merger, or efficient radial migration of stars. 
Subject headings: Galaxy: abundances — Galaxy: disk — stars: abundances — stars: late-type 



1. INTRODUCTION 

Since its identification (|Gilmore fc Reidl 11983ft . the 
thick disk has been shown to have distinct kine- 
matics (e.g. IChiba fc Beers! [20001 iGilmore et all [20021: 
iSoubiran et alj 12003ft aiula distinct metallicity distri- 
bution (e.g. iMaiewskil 119931: IChiba fe Beers! [2000) from 
the stellar halo and thin disk of the Milky Way Galaxy. 
The formation history of the thick disk can provide con- 
straints on the origins and formation of the thin disk 
and halo, and ultimately the Milky Way Galaxy itself. 
The formation mechanism of the thick disk, however, has 
been the subject of much discussion for decades. 

Soon after the discovery of the thick disk, 
Uones fc Wvsd (|1983l ) proposed that early star for- 
mation in a thin disk, which formed before the Galactic 
potential reached virial equilibrium, could be subse- 
quently flattened by violent relaxation to form the 



thick disk. Later, iBurkert et al.l (|1992ft showed that 
a thick disk could be formed as the result of a rapid, 
dissipational collapse accompanied by star formation. 

The first model to gain traction, however, was that 
in which a pre-existing thin stellar disk is heated by 
a merger with a fairly robust and mass ive satellite, 
around ~ 20% the mass of th e disk (e.g., lOuinn et al.1 
[l993HVelazauez fc Whiti[l999[ ). The thick disk is then 
formed primarily from the stars in the 'thickened' disk, 
while some stars are directly accreted from the satellite 
galaxy. More recently, this model was extended to in- 
clude cosmologically motivated heating by accretion and 
merging from many sa tellites (lHavashi fe Chibal 120061: 
iKazantzidis et alj 120081: iVillalobos fc Helmi! 12008ft . how- 
ever, the most massive satellites still dominate the heat- 
ing. 

Another model has the thick disk forming from the di- 
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rect a ccretion of stars from satellite galaxies ([Abadi et al.1 
2003). Alternatively, the thick disk could have formed 
during a period of rapid star format ion associated with 
multi ple minor, gas-rich mergers (Brook e t al.l 12004 
2005) . In this model, the majority of the thick disk stars 
are formed from the accreted gas in the primordial disk, 
with a small minority being directly accreted during the 
mergers. 

Models that have more recently gained traction are 
those that do not necessarily require a merger history 
to the formation of the thick disk. Instead, the thick 
disk could have formed from stars that have radially mi- 
grated outward from the inne r disk by resonant scatter- 
ing by transient spiral arm s (ISellwood fc Binnevl [2002; 
iSchonrich fc Binnev 2009a. b). Another model, proposed 
bv iBournaud et al.l <j2009h has the thick disk forming 
from internal gas clumps in the disk at high redshift. 

What is the best way to distinguish these models? The 
distributions of kinematics and metallicity provide some 
discriminants, and these differ between the models most 
strongly at the metal-poor end. While there is clearly not 
a simple one-to-one relationship between age and metal- 
licity, the most metal-poor stars in general formed early, 
in every scenario. 

Elemental abundance patterns reveal even more de- 
tailed information about the star formation history and 
chemical evolution of a stellar population. The ratio of 
a-elemental (e.g., Mg, Si, Ca, and Ti) abundances to 
iron abundance is sensitive to the relative numbers of 
core-collapse Type II supernovae (SNe II) and Type la 
supernovae (SNe la) that have occurred in the past. The 
a-elements are primarily synthesized in massive stars and 
ejected in SNe II on short timescales (~ 10 7 yr) after the 
formation of the progenitor star, with the ratio of [a/FeJ3 
from any one SNe II depending on the mass of the pro- 
genitor massive star (|Kobavashi et all 12006ft . SNe la, 
which result from white dwarf remnants in intermedi- 
ate mass binary systems with mass transfer, are major 
sources of iron (and not a-elements) and occur on longer 
time scales than SNe II. The actual distribution of explo- 
sion timescales depends on the details of the model for 
the progenitors of Type la SNe, but the onset is always 
later than that of Type II SNe, and there are always ex- 
plosions several Gyr after the initial star formation (e.g. 
Matteucci et al. 2009). 

A self-enriching system will then show enhanced [a/Fe] 
ratios (compared to the Sun) early in the star-formation 
process, when enrichment is dominated by core-collapse 
SNe II, and the actual enhancement is determined by 
the initial mass function (IMF) of massive stars (assum- 
ing good mixing of the interstellar medium, ISM, and a 
well-sampled IMF) . Stars formed after significant contri- 
butions from SNe la to the ISM (~ 10 8 - 10 9 yr after the 
first star formation episode) will have considerably more 
iron enrich ment, and so they will have decreased [ct/Fe] 
ratio s (e.g.JMatteucci fc Bro cato 1990t lGilmore fc Wvsd 
119911: lMatteuccill2003h . ~ 

The models can therefore be tested using a sample of 
metal-poor ([Fe/H] < —1) thick disk stars that probe a 
volume out to several kiloparsec around the Sun. High 
resolution spectroscopic observations of each star is nec- 

1 We use bracket notation to indicate relative abundance ratios of 
two elements A and B: [A/B] = log[n(A)/n(B)] -log[n(A)/n(B)] . 



essary to determine the detailed abundances of the a- 
elements in each star. We identified for study a sample 
of candidate metal-poor stars in the thick disk selected 
from the Radial Velo city Experiment survey (RAVE, 
iSteinmetz et al.1 120061 ) , and observed this sample with 
high resolution echelle spectrographs at several facilities 
arou nd the world. 

In lRuchti et al.l (|2010t hereafter R10), we derived sev- 
eral [a/Fe] ratios for a subsample of red giant branch 
(RGB) and red clump/horizontal branch (RC/HB) stars 
from our full candidate metal-poor thick-disk sample. 
We found a significant fraction of these stars (~ 40%) 
to be most consistent with the thick disk. The abun- 
dances of these metal-poor thick-disk stars had high 
[a/Fe] ratios, consistent with rapid star formation, 
which agrees with previous studies o f much smaller 
samples of metal-poor thick-disk stars (Fulbrieht 2002; 
i Bensbv et al.ll2003t iBrewer fc Carnevll2006l: iReddy et al.1 
l2006t iReddv fc Lambert! 120081: lAlves-Brito et al.M2010D . 
Further, we found that the [a/Fe] ratios were indistin- 
guishable from those of the halo, indicating that the halo 
and the thick disk shared a similar massive star IMF and 
similar efficient mixing of enriched material into the in- 
terstellar medium. 

The high [a/Fe] ratios of the metal-poor thick-disk gi- 
ants provided constraints on those models of the forma- 
tion of the thick disk that are driven by direct accre- 
tion of stars from satellite galaxies. Stars in present-day 
Milky Way satellites have lower values of [a/Fe] at a 
given [Fe/H] (for [Fe/H] > -2) than th e stars of the 
stella r halo, thick disk, or thin disk (see iTolstov et al.1 
2009) . This is understood in terms of the differences in 
the star formation histories of the stellar populations (as 
described above), with the dwarf galaxies having very 
slow enrichment, and strongly non-monotonic star for- 
mation rates. In the models that include direct accretion 
of stars, accretion takes place well beyond 1 Gyr after 
the initial star formation episode. This implies that if 
each satellite had star formation rates similar to surviv- 
ing dwarfs, then the accreted stars would have formed 
after significant iron contribution from SNe la, and thus 
would have low [a/Fe] ratios. From this we concluded 
that dwarf galaxies similar to present-day dwarf galaxies 
did not play a major role in the formation of the thick 
disk. 

Further analysis, however, is needed to distinguish 
those models that do not include a significant fraction 
of accreted stars from satellite galaxies. These models 
all have early, rapid star formation, consistent with high 
[a/Fe] ratios. A useful diagnostic is the amplitude of the 
radial and vertical abundance gradients predicted by the 
models. For example, the gas-rich merger model predicts 
more uniform chemical abundances, while low-amplitude 
vertical gradients are possible in the heating scenario. 
Slow, dissipational settling, on the other hand, would 
produce a significant vertical gradient in metallicity, as 
well as a possible vertical gradient in the [a/Fe] ratios. 

In this work, we extend the analysis from R10 to our 
entire sample of metal-poor thick-disk candidates. We 
further augment that study with a detailed analysis of 
the abundance gradients in the thick disk. In §2 and §3, 
we briefly describe the candidate selection and high res- 
olution spectroscopic observations. In §4 we describe our 
stellar parameter analysis. We give a full description of 
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our distance estimation techniques in §5, and briefly dis- 
cuss our population assignments in §6 (a full description 
can be found in RIO). In §7 we report on the abundance 
correlations and gradients seen in the data for our full 
sample, and show that our conclusions from RIO have 
not changed by including the dwarfs and sub-giants in 
our sample. We further quantify our abundance results 
using IMF-weighted yields in §9. We then use orbital 
eccentricities of our sample stars in §10 to distinguish 
further the models of the formation of the thick disk. 
Finally, we discuss our findings in §11 and conclude in 
§12. 

2. SELECTION OF CANDIDATE STARS FOR STUDY 

As in RIO, all candidate stars were first selected from 
the internal RAVE catalog (Version DR2) to have cali- 



M/Hl ca i < —0-7 (fo r more details on 
Zwittcr ct al. 200§) and T eS values 



brated metallicities 
the calibration, see 
between 4000 K and 6500 K. The low temperature cut 
was applied to reduce contamination from metal-rich gi- 
ant stars and M-dwarfs. Stars hotter than 6500 K have 
a higher likelihood of being rapid rotators and may also 
have larger non-LTE effects, which could affect our ele- 
mental abundance analysis (as described below). Those 
stars that met these parameter constraints and had the 
highest probability of being thick disk stars according to 
their kinematics (based on the RAVE parameters and our 
distance estimates, for which the procedure is described 
below) were selected for high resolution spectroscopic ob- 
servations. 

3. HIGH-RESOLUTION ECHELLE OBSERVATIONS 

Observations were conducted during the period be- 
tween 2007 May and 2009 February. All high resolu- 
tion spectra were obtained u sing one of the follow ing 
echelle spectrographs: MIKE ([Bernstein et al.l 120031 ) on 
the Magellan- Clay tele scope at Las Camp anas Observa- 
tory in Chile, FEROS (JKaufer et al.lU999h on the MPG 
2.2-m te lescope at ESQ La Si lla Observatory in Chile, 
UCLES ([Walker fc Diegoll 19851) on the Anglo- Australian 
telescope in Austra lia, and the ARC echelle spectrograph 
(jWang et al.ll2003l ) on the Apache Point 3.5-m telescope 
in New Mexico, this last facility for stars visible from the 
northern hemisphere. 

The instumental set-ups gave resolving powers between 
35,000-45,000, providing complete spectral coverage from 
A3500A to A9500A, except for UCLES which had com- 
plete spectral coverage between A4460A and A7270A. 
The echelle spectral data were reduced following the 
same procedures as described in R10. The final reduced 
spectra all yielded an S/N ratio greater than 100 pixel -1 
at A5000-6000A and a minimum S/N - 40 around 
A4000 A, which is sufficient for detailed elemental abun- 
dance analysis. 

3.1. The Final Sample 

Our full sample of metal-poor thick-disk candidate 
stars consists of the sample from R10 (212 RGB stars 
and 31 RC/HB stars) and 74 main-sequence or sub-giant 
(MS/SG) stars. An additional 2 RGB stars were discov- 
ered during the analysis of the MS/SG stars, bringing our 
total to 319 candidate metal-poor thick disk stars with 
high resolution spectroscopic observations. As indicated 
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Figure 1. Distribution of the difference between the echelle- 
derived heliocentric radial velocities and those estimated by RAVE. 
The average difference is — 1 ± 3 km s — 1 . 

in R10, ten of these candidates (5 RGB and 5 MS/SG) 
were observed twice for internal consistency checks. Ta- 
bled] lists the observational data for the sample. 

The full sample has heliocentric radial velocities de- 
rived from the echelle spectra that differ from those in 
the RAVE database by an average — 1 ± 3 km s -1 . Fig- 
ure Q] shows the distribution of the differences in the ra- 
dial velocity estimates. The standard deviation of the 
RAVE radial velocities is of the order of 2.5 km s -1 , 
which is slightly lower than that of our comparison. 
We cross-checked our final sample with the list of spec- 
troscopic binar y stars found in the RAVE survey by 
IMatiievic et al.l (|2010D . and found no matches. Those 
stars with differences of 10-15 km s _1 could possibly be 
in low-amplitude, single-lined binary systems in which 
the primary to secondary ratio is about 5-10 to one. This 
type of binary system would not have a significant effect 
on our subsequent analysis and results and we retained 
these stars. 

The coverage in Galactic (l,b) coordinates of our final 
sample is illustrated in Figure [2) RAVE mainly targeted 
fields at Galactic latitudes greater than 20° (only target- 
ing a few low-latitude fields), which is why there are very 
few stars with \b\ < 20° in our sample. 

It is clear that our sample of metal-poor thick-disk can- 
didates does not have the same completeness as the entire 
RAVE sample from which we selected our candidates. 
Figure [3] compares the distribution of the RAVE internal 
/-magnitude of our sample to that of the entire RAVE 
sample. Most of our stars lie between 9 < I < 11, about 
a half magnitude brighter than typical for the parent 
RAVE sample. In addition, our selection function was 
by no means homogeneous during the observing runs. It 
is important to realize, then, that our sample does not 
satisfy either magnitude- or volume-completeness. 



4. ABUNDANCE ANALYSIS 

Initial estimates of the stellar parameters (T- ln - u log 
[Fe/H]i n i) for each star were determined foll owing the 
methodology of R10, based on the methods of Fulbright 
(|2000. hereafter FOOl ). Full details of this technique can 
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Table 1 

Observational Data 



Star 


RA° 


DEC 




Obsdate 6 ^ 


Observatory 


S/N d 




(°) 


(°) 




(yyyymmdd) 




(pix" 1 ) 


C0023306-163143 


5.878 


-16.529 


11.2 


20081015 


LCO 


160 


C0315358-094743 


48.899 


-9.796 


9.9 


20081016 


LCO 


170 


C0408404-462531 


62.169 


-46.425 


11.9 


20081016 


LCO 


190 


C0549576-334007 


87.490 


-33.669 


11.0 


20081015 


LCO 


170 


C1141088-453528 


175.287 


-45.591 


10.3 


20070506 


LCO 


100 


abed 





























Note. — Tablc^is published in its entirety in the electronic edition of the Astrophysical 
Journal. A portion is shown here for guidance regarding its form and content. 

a equinox 2000 

b 2009combo or 2007combo refers to the combined spectra from two different nights of 
the same run. 

c For APO runs, 'combo' refers to combined spectra from different runs. 
d Estimated between A5500-6000 A. 
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Figure 2. Sample coverage in Galactic (l,b) coordinates. The 
dotted line represents the Galactic plane, while the dashed curve 
illustrates the position of the celestial equator in Galactic coordi- 
nates. 
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Figure 3. Distribution of /-magnitude for our candidate sam- 
ple (red-dashed) and the parent RAVE sample (solid-black). Note 
that the RAVE sample only includes those stars that satisfied our 
constraints in [[2] 

be found in RIO. It is important to note, however, 
that this method computes the initial quantities using 
one-dimensional, plane-parallel LTE Kurucz model at- 
mospherefl 

2 http://kurucz.harvard.edu/ 



It is important that we derive accurate stellar param- 
eters since we later use these parameters to determine 
distances to our stars. Parameters derived by our initial 
analysis, however, can often show large differences from 
expect ed values compute d by different methods. For ex- 
ample, llvans et al.l (|2001h show, using giants analyzed in 
the relatively metal-poor ([Fe/H] ~ —1-2) globular clus- 
ter M5, that the log g derived by methods similar to our 
initial analysis is often too low, especially for stars near 
the RGB-tip. This is in part due to the fact that Fei 
(and somewhat Fen) can be strongly affected by non- 
LTE effects, and so the estimate of surface gravity and 
iron abundance can be quite unreliable, especially in low- 
gravity giants and metal-poor stars when ionization equi- 
librium b etween Fei and Fell is assumed in our initial 
analysis dThevenin fc IdiartJ 119991: lAsplund et all 119991: 
Mashon kina et al. 1120111 Lind, K. & Ber gemann, M., pri- 
vate communication). As summarized in RIO, we there- 
fore utilized high-resolution spectroscopic data for sev- 
eral globular cluster stars, as well as reanalyzed several 
F00 RGB stars with good Hipparcos parallaxes, to test 
the accuracy of our initial stellar parameter analysis (de- 
rived from our echelle-based analysis) for the giant stars 
(RGB and RC/HB) in our sample. Using independent 
estimates of the stellar parameters for these test cases, 
we found that the estimates from our initial technique 
must be systematically corrected to achieve improved ac- 
curacy. Using a sample of 28 F00 MS/SG stars, we have 
also found that the initial parameters for our MS/SG 
stars must be similarly corrected. In the following sec- 
tions, we give a full description of the corrections to the 
stellar parameters for our entire sample. 

4.1. Giant Star Parameters 

In R10, we provided a brief summary of the correc- 
tions needed for the giant stars in our sample. In this 
section, we give the full description of those corrections. 
It is important to note that the giant star test samples 
did not include RC/HB stars, however, since the stellar 
parameter corrections are independent of the evolution- 
ary state of the giant star, it was assumed in RIO that 
the RC/HB stars would follow the same corrections as 
the RGB stars. 

We c omputed a photometric temperature (T P hot), us- 
ing the lGonzalez Hernandez fc Bonifaciol (|2009f) 2MASS 
color-temperature transformations, and a 'bolometric' 
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gravity ((jw = AttGMctT 4 / L) for each globular cluster 
and F00 RGB star. A photometric temperature was also 
computed for the stars in our sample that showed low 
initial reddening (Aj < 0.05, see §5.4|) . We found in RIO 
that the offset, At = T ini — T p i lot , shows different cor- 
relations with [Fe/H]i n i depending on the initial Tj n j of 
the star. The behavior can be well approximated by two 
regimes, hotter and cooler than 4500 K, with correlations 
as illustrated in Figures [4] and [5j 

A robust least-squares linear fit to the stars with Tj n i > 
4500 K (see Figure 0J resulted in: 







A 



if [Fe/H] ini > -1.0 
240 if - 1 > [Fe/H] ini > -2.5 
if [Fe/H] ini < -2.5 

(1) 

while for stars with T ini < 4500 K (see Figure [5]): 



240 [Fc/H] 
-400 



A T = 



if [Fe/H]. 

120 [Fe/H] iui + 140 if [Fe/H] 
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Figure 4. At vs. Ti n i, log gi n i, and [Fe/H]i n i for giant stars with 
Tini > 4500 K. Black circles and red triangles represent the globular 
cluster stars and F00 RGB stars, respectively. The blue x-symbols 
are giants from our sample that have small initial reddening (Aj < 
0.05). Notice that there is no obvious trend with Ti n i or log cji n i, 
however, it is clear that T"i n i — T pnot depends on [Fe/H]i n i. 

Following the reasoning of R10, we adopted the above 
corrections, adding At to T; n i to obtain final values for 
the stellar temperature (denoted as T e g), to bring our 
initial temperatures to a scale that will reduce spuri- 
ous trends as found for our initial spectroscopic analy- 
sis. We then used this final estimate of T e g in the abun- 
dance analysis to obtain a new (ionization-balanced) es- 
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Figure 5. Same as Figure[4] but only for stars with Ti n \ < 4500 K. 
Similarly to Figure[4] the only obvious correlation is with [Fe/H]i n i. 
It is important to notice, however, that the slope of the correlation 
is shallower than for the stars with T; n ; > 4500 K. This is the 
reason we separate the stars into two temperature regimes. 

timate of gravity, denoted log <7 p hot- Comparisons of 
this log 5 P hot estimate with log gboi showed improve- 
ment, but an offset still remained for the lowest gravity 
stars. This residual offset correlates with log <7 p hot (see 
Figured]) such that a least-squares linear fit to the data 
resulted in: 



A 9 = 



0.0 

0.6 log .gphot - 0.6 



if log # p hot > 1.0 r „N. 
if log «? phot < 1.0 [6) 

For those stars with log g p hot < 1, we applied the 
above correction to the log g p hot estimate. We then 
adopted both the corrected T e g and this new log g esti- 
mate to get a new estimate of the iron abundance. We 
chose the iron abundance from Fe II as our final esti- 
mate of iron abundance to reduce sensitivity to non-LTE 
effects (jThevenin fe Idiartl [l99l lAsplund et all fl999h . 
Scatter from our final estimates of temperature and grav- 
ity around T p hot and log <?boi provided error estimates, 
OT eff = 140 K and a\ og g = 0.2 dex. The iron abun- 
dances in the literature are not on a uniform scale, so we 
estimated the error crpo/H] =0.1 dex, from star-to-star 
scatter within any one globular cluster. 

4.2. MS/SG Star Parameters 

Only the giant stars in our sample were tested in R10. 
We therefore reanalyzed 28 MS/SG stars from F00 to 
test the accuracy of the initial stellar parameters for the 
MS/SG stars in our sample. We computed an indepen- 
dent log <7boi and photometric temperature, Tphot, fol- 
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Figure 6. Difference, log g p hot — log 9bol, vs - corrected T e fj, 
log g p hoti an< i [Fe/H]ini. Black circles and red triangles represent 
the globular cluster stars and F00 RGB stars, respectively. Notice 
that after the temperature correction (equations [l] and [2]| A 9 only 
depends on log g pho f 

lowing the same methods as for the giants. We also 
included estimates of photometric temperature for the 
MS/SG stars in our sample with low initial reddening 
(as done for the giant stars). We found, through a linear 
fit to the data, that the offset, At = Tlni-Tphot, correlates 
with [Fe/H]i„i for these stars such that: 



A T = 190 [Fe/H] ini + 150 



(4) 



This correlation is illustrated in Figure [JJ In addi- 
tion, Figure [5] shows that the difference, A g = log — 
log <7bob correlates with both [Fe/H]i„; and log for 
the MS/SG stars. The linear fit to this correlation is 
given by: 



A g = 0.3 [Fe/H] ini + 0.2 log g ini - 0.8. 



(5) 



We simultaneously corrected our initial T ml and log 
values according to the above correlations, which resulted 
in mean scatter around T p h t and log (/boi of 150 K and 
0.2 dex, respectively. Again (as for the giants), we chose 
the abundance of Fe n as our final estimate of iron abun- 
dance. 

4.3. Final Values of Stellar Parameters 

We used the above procedures to ensure that the 
derived parameter values for our candidate metal-poor 
thick disk sample are on the same scales as the globular 
cluster stars and Hipparcos stars. Stars with repeat ob- 
servations showed small differences in the values of the 
stellar parameters after the corrections above, but they 




Figure 7. Difference, T lni — T phot , vs. T ini , log g ini , and [Fe/H] ini 
for the 28 MS/SG stars from F00. Red triangles represent the 
F00 data, while the blue x-symbols represent MS/SG stars in our 
sample with low initial reddening (Aj < 0.05). 



are all smaller than our estimated errors during the cor- 
rection procedure. The effective temperatures showed 
mean differences of 16 ± 58 K, the log g values differed 
by 0.05 ± 0.12 dex, and there was a mean difference of 
0.02 ±0.07 dex in [Fe/H]. 

Table [2] gives our final stellar parameter values for our 
sample. Figure [9] illustrates the position of our stars on 
the log g versus log(T e ff) plane before (left panel) and af- 
ter (right panel) the above correcti ons. Two isochrones , 
computed by the Padova group (jGirardi et all 12001 
Marigo et al. 2008), with an age of 12 Gyr and metallic- 
ities Z — 0.001 and Z — 0.006 are plotted for reference. 
Notice that the resulting gravities for the RGB trace the 
isochrones much better after the corrections. Note that 
some MS/TO stars have different metallicities than the 
reference isochrones, which makes them appear to be in 
an impossible part of the plane. 

In Figure 1101 we compare our stellar parameter val- 
ues with those from RAVE. As can be seen in the figure, 
on average the RAVE values for T e g and log g match 
the echelle-derived values. There is, however, a large 
spread, especially for RGB stars at low (log g < 2) grav- 
ity. The echelle-derived values also tend to have hotter 
temperatures and larger gravities for the dwarfs and sub- 
giants. Metallicity comparisons show a tight trend be- 
tween the difference of the two measurements with the 
echelle-derived value. It is clear that the RAVE [M/H] 
is not the same as our echelle-derived [Fe/H]. 

4.4. Abundances of the Alpha- Elements 
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Figure 8. A 9 vs. corrected T e g, log g; n ;, and [Fe/H]; n ; for the 
F00 MS/SG stars. It is clear from this figure that A 3 for the 
MS/SG stars depends on both log g- m \ and [Fe/H]i n j. 
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Figure 9. Stellar positions in the log g vs. log(T e ff ) plane before 
(left panel) and after (right panel) our stellar parameter correc- 
tions. The isochrones plotted have an age of 12 Gyr and metallici- 
ties Z = 0.001 (solid red curve) and Z = 0.006 (dashed blue curve). 
Some MS/TO stars have different metallicities than the reference 
isochrones, which make s th em appear to be in an impossible part 
of the plane. See Figure [T2l for plots of isochrones more appropriate 
for the MS/SG stars. 



Elemental abundances of several a-elements were 
taken from R10 for the RGB and RC/HB stars, while 
those for the MS/S G stars were de rived using the MOOG 
analysis program (Sn edenl 119731 ) using the final stel- 
lar parameters for each star (as done for the giants in 
R10). The ratio of the abundance of each a-element to 
iron abundance , comp uted using solar abundances from 
iGrevesse et al.l ()1996[ ). can be found in Tabled! Com- 
parisons of repeat observations for each ratio, showed 
differences of ^ 0.03 ±0.1 dex. We therefore adopted 






Figure 10. Stellar parameter comparisons with RAVE values: 
T e ff (upper left panel), log g (upper right panel), and metallicity 
(lower left panel). AH panels plot the echelle value of the pa- 
rameter minus the analogous RAVE parameter versus the echelle 
value. Note that there is a tight correlation between the differ- 
ence in metallicity values versus the echelle-derived [Fe/H], which 
results from the fact that the RAVE [M/H] is the estimated total 
metallicity from iron and the alpha elements together. The more 
metal-poor stars are more enhanced in the alpha-elements (shown 
later), which increases this total metallicity. 

an error of 0.1 dex in the [a/Fe] ratios. Final abundance 
values for stars with repeats were taken as the average 
of the two estimates for each element. 



5. STELLAR DISTANCES 

For this effort to yield meaningful conclusions, we have 
to determine meaningful population assignments (see 
Sj6]). Meaningful population assignments require accu- 
rate kinematic values to distinguish the thick disk from 
the halo and thin disk. The greatest contributor to the 
uncertainty in the kinematics is the distance uncertainty. 
We therefore applied techniques that minimized the un- 
certainty in the derived distances to our stars. 

5.1. RGB Stars 

Distances to the RGB stars were adopted from R10, in 
which we utilized a method that determines probability 
weights for a set of isochrones. An absolute Mk s mag- 
nitude of a given star was then determined by matching 
the parameters of a star to the most probable isochrone 
using the technique described below. 

5.1.1. The Isochrones 

Throughout the analysis of R10, as well as this work, 
we used the set of isochrones computed by the Padova 
group. There are, however, many other isochrone sets 
available. To test the systematics resulting from using 
one set of isochrones over another, we als o utilized a sec- 
ond set of isochrones , the Yonsei-Yale ()Yi et al.l 120011 : 
iDemarque et al.l I2004D isochrones, in the analysis of a 
small subset of our stars covering the full range of stellar 
parameters for our entire sample. The distance estimates 
derived from fitting the test data to each isochrone set 
resulted in an average difference in distance of about 10% 
- well inside our experimental error (as will be discussed 
below). This implies that isochrone differences will not 
have a noticeable effect on our estimated distances. 
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Table 2 

Atmospheric and Abundance-Ratio Data 



Star T cff (K) log g [Fe/H] a v t Mg 6 Si 6 Ca 6 Til 6 Till c 



C0023306-163143 


5528 


3.29 


-2.30 


1.4 


0.32 


0.46 


0.34 


0.11 


0.15 


C0315358-094743 


4722 


1.77 


-1.35 


1.6 


0.30 


0.45 


0.28 


0.11 


0.24 


C0408404-462531 


4600 


0.86 


-2.13 


2.1 


0.42 


0.44 


0.32 


0.09 


0.18 


C0549576-334007 


5379 


3.04 


-1.76 


1.3 


0.24 


0.22 


0.31 


0.15 


0.21 


C1141088-453528 


4592 


1.01 


-2.32 


1.9 


0.35 


0.37 


0.28 


0.02 


0.20 



Note. — Tabic [2] is published in its entirety in the electronic edition of the Astrophysical 
Journal. A portion is shown here for guidance regarding its form and content. 

a given as [Fell/H] 
b given as [X/Fel] 
c given as [X/Fell] 



The Padova isochrones were chosen for the remain- 
der of our analysis since they can be simply computed 
for 2MASS magnitudes, and also include the late stages 
(HB and AGB) of stellar evolution. For the RGB dis- 
tance analysis, we created a grid (without interpolation) 
of Padova isochrones with metallicities ranging between 
Z = 0.0001 and Z = 0.01 (with solar elemental ratios) 
with a step-size of 0.0002, and ages ranging between 
4—12 Gyr with a step-size of 1 Gyr. 

5.1.2. The Fit Technique 

The technique from R10 fit the stellar parameters of 
each star to several isochrones. A series of sequential 
weighted averages over the isochrones were then per- 
formed to find the most probable Mk s magnitude for 
each star. Note that we do not perform any interpo- 
lation between isochrone grid points. The points along 
the RGB of each isochrone are approximately uniformly 
distributed, however, the spacing between points on dif- 
ferent isochrones (especially of differing metallicity) can 
be uneven, which may cause over-sampling or under- 
sampling of different regions of the entire grid. We at- 
tempted to reduce this effect by restricting the range of 
metallicities used in the fit and introducing a prior prob- 
ability to each grid point, as described below. 

The parameters used to fit each star were T e ff, log g, 
and the Z-metallicity of the star. The Z-metallicity was 
determined by combinin g |Fe/H] and the m ean [a/Fe] 
using the prescription of Sa laris et al.l ()1993t ). 



Z = Z o (0.638 • 10 



[a/Fe] 



0.362), 



(6) 



where Z a = Z o 10l Fe/H1 and Z Q = 0.019 is the solar 
metallicity. Errors in [Fe/H] and [a/Fe] were propagated 
to estimate the errors in Z. Next, the probability distri- 
bution of each parameter was assumed to be described 
by a Gaussian, 



Af(X*,a x J oc exp[-(X - X+f /2a 2 x ] 



(7) 

The sigmas, ox„ i are the estimated errors in each pa- 
rameter, X+ — {T c ff, log g, Z}, derived for each star in 
fjU The grid of isochrones were then redefined so that 
for a single age, the grid is limited to only isochrones 
with metallicities within ±10 times the error of a star's 
Z-metallicity. 

In an effort to reduce cases in which two isochrones of 
differing ages and metallicities might give the same prob- 
ability, a prior probability function was also introduced. 



This function, computed for each point in the grid of 
isochrones, was derived from the BaS TI luminosity func- 
tion tracks (jPietrinferni et ahl [2004) . By this method, 
evolutionary stages with longer lifetimes will be assigned 
a larger prior probability than short-lived stages, thus 
biasin g the fit to the longer-lived stages (cf. iPont fc Everl 
l200l . 

The probability information was then combined so that 
for each point on an isochrone of a given age, i, and 
metallicity, j, a probability weight was computed by, 



(8) 



where J\f(X+, ax„ ) is given by equation ([7]), ^^(Xk) is 
the prior probability, and X k equals the set of values, 
{T e ff, log g, Z}, at the given point, k, on the isochrone. 

The total probability weight for an isochrone of a 
given age and metallicity was determined by summing all 
Pij(Xk) over all k points in the isochrone. Further, the 
most probable Mr- s for the given star when matched to 
that isochrone was computed by performing the weighted 
average, 

E k Pv(Xk)M l3 (X k ) 



Mi. 



(9) 



where My is the most probable absolute Mk s magnitude 
of the star if it has the age and metallicity of the given 
isochrone. 

Next, the weighted average of the estimates, My, of 
absolute magnitude was performed over all j-metallicities 
to give the most probable estimate Mj of a given age, 



Mi 



T.3P* 



(10) 



where Pj = J^j Pij 1S the total probability weight at i- 
age, summed over all metallicities within the range. Fi- 
nally, the estimates of the best absolute magnitude, Mj, 
obtained for isochrone ages of 10, 11, and 12 Gyr, were 
combined by a weighted average to produce the final es- 
timate of the absolute magnitude Mk s for each star (as 
used in R10). 

5.2. RC/HB Stars 

As with the RGB stars, distances to the RC/HB stars 
were adopted from R10. Mass loss along the RGB, which 
affects the position of the HB, is not well modeled in the 
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Padova isochrones. The absolute Mk s magnitude for 
each RC/HB star was therefore assumed to be equal to 
the single HB point on the isochrone of equivalent Z- 
metallicity and age = 12 Gyr. 

5.3. MS/SG Stars 

Distances determined by fitting near the turn-off (TO) 
on isochrones are much more sensitive to age determina- 
tion than on the RGB. In addition, the fit technique for 
the RGB stars described above can result in unphysical 
solutions in this region. The percent error in the T e g 
value of a star is typically much smaller than that of the 
log g values, implying the probability distribution close 
to the TO region can be double-peaked. For example, a 
star with a gravity that lies to the right of the TO point 
on an isochrone will have a probability peak on the MS- 
branch as well as one on the SG-branch. This will cause 
any weighted averaging of the isochrones to choose the 
most probable fit to be a point not on the isochrone itself, 
which is physically impossible. We, therefore, employed 
a different method for finding distances to our MS/SG 
stars to reduce this effect. 

Instead of fitting to isochrones, we estimated the ab- 
solute Mks magnitude by assuming a mass, using the 
definition of effective temperature and adopting a bolo- 
metric correction for each MS/SG star, which is given 
by: 

/ M 
M K „ =-2.5 log 

+4 log - log -L ) + A/boi,o - BC Ks 
Tq 9qJ 

where T = 5770 K, log g G = 4.44, M bo i, = 4.72 is the 
absolute bolometric magnitude of the Sun, and BCk are 
the b olometric corrections derived from iMasana et al.1 
(2006), which are applicable for the entire metallicity 
range of our stars and have propagated errors around 
~ 0.05. This equation, however, still depends on the 
mass of a star. Investigating the mass ranges around the 
turn-off on the isochrones shows that MS/SG stars range 
in mass between ~ 0.8 — 1.0 M© for ages 4 — 12 Gyr. 
Adopting a mass of 0.8 M implies that the star is old, 
while adopting a mass of 1.0 M implies that the star is 
younger. Comparing a star's position in the log g ver- 
sus T e g plane to a 4 Gyr and 12 Gyr isochrone, we can 
determine a reasonable mass estimate. This is described 
in more detail below. 

5.4. Reddening 

An estimate of the reddening to a given star was 
computed during the proced ure to obtain a distan ce by 
the following method. The iSchlegel et all (|1998f ) dust 
maps and extinction calculator were initially used to 
calculate reddening to infinity in a specific line-q f -sight . 
iBonifacio et~aH (|2000f ) found that ISchlegel et all (|1998f ) 
overestimate the reddening when E(B — V)s c h > 0.10 
, and 52 of our candidate stars meet this criterion. We 
therefore adopted their correction for the extinction: 

E{B~V) C = 

f E(B - V) Sch , if E(B - V) Sch < 0.10 
\0.10 + 0.65 [E(B - V) S ch - 0.10], otherwise. 

(12) 



where E(B — V)c is the reddening estimate after the 
correction. Several stars (especially MS/SG stars) are 
close enough to lie within the Galactic dust distribution, 
and thus the reddening should be reduced. Adopting 
a simple exponential model for the dust, with a scale- 
height of h — 125 pc (jBonifacio et al.lfeOOOf ). the redden- 
ing to a star at distance D and Galactic longitude b is re- 
duced by a factor 1 — exp(— \D sin b\/h). The reddening 
was recomputed iteratively until the difference between 
the current estimate of t he distance and the p revious one 
was less than 2 percent. Schlcgcl et al.l (|19 98) assumed a 
Rv = 3.1 extinction curve to compute extinction coeffi- 
cients in other passbands. They do not, however, derive 
values for the 2MASS passbands. We determined our 
extinction coefficients by assuming the extinction coeffi- 
cients fromlMcCajJ (pOOl . 

Aj = 0.819 E(B - V) c ri „v 
A Ks = 0.350 £(B - T/) c . 1 j 

The mean values of Aj and Ak s for our sample are both 
less than 0.1, with errors less than 0.1 mag. An error 
of 0.1 mag in reddening will only cause a ~ 4% shift in 
distance, which is within our error estimates described 
below. 

5.5. Comparisons & Error Analysis 

We adopted the estimate of 20% error on the distance 
to RGB and RC/HB stars from R10, which was based 
on comparisons between our estimates and literature esti- 
mates for the globular cluster stars and estimates for the 
F00 RGB stars derived from Hipparcos parallaxes. For 
the MS/SG stars, we applied our technique from above 
to estimate distances for each of the 28 F00 MS/SG stars, 
and compared to the Hipparcos distances. Figurc[TT]plots 
these stars in the log g versus log(T e ff) plane. The ma- 
jority of the stars appear to be old. There are, however, 
several stars with [Fe/H] > —0.6 that could have younger 
ages. If we assume all of the stars are old, and choose 
a mass of 0.8 M , then the distance estimates from our 
technique differ from those derived from the Hipparcos 
parallaxes by only 1% ± 17%. Increasing the mass to 
0.9 M and 1.0 M increases the difference to 7%± 18% 
and 13% ± 19%, respectively. It is important to notice 
that if we increase the mass to one that is intermediate 
between older and younger ages, then the difference only 
increases by 6%, which is well within the scatter. 

Figure [T2l shows the same plots as Figure mi except we 
now plot our metal-poor sample MS/SG stars. It is clear 
that all stars with [Fe/H] < — 1 are old. We therefore as- 
sume they have a mass of 0.8 M . We assume a mass of 
1.0 M for the two stars with [Fe/H] = (-1, -0.8) close 
to the 4 Gyr isochrone (log(T e ff) > 3.78), while we adopt 
0.8 M for the remaining stars in that bin. Those stars 
with [Fe/H] > —0.8 show no clear separation in age. We 
therefore chose an intermediate mass of 0.9 M for all 
MS/SG stars with [Fe/H] > -0.8. If we apply the above 
prescriptions to the F00 MS/SG stars, we find a differ- 
ence between our bolometric distances and the distances 
derived from parallaxes to be only 4% ± 17%. We there- 
fore adopt a 20% error on the distances to our MS/SG 
stars (similarly to the giants and RC /HB stars) , so as to 
include both offsets and scatter. 

Distance estimates based on RAVE pi peline values of 
stellar parameters are also now available ( Brcddels et al. 
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Figure 11. Plot of log g vs. log(T cff ) for the F00 MS/SG stars 
(after all parameter corrections) for different metallicity bins as 
shown in each panel. In addition two isochrone curves are plotted 
(for the given Z-metallicity shown in each panel). The black curve 
and red curve represent ages of 12 Gyr and 4 Gyr, respectively. 
Note that the majority of the F00 stars appear to be old, except 
for a few with [Fc/H] > -0.6. 

[20101: iZwitter et all [20Toi: iBurnett fc Binnevl [20loh . 
Our distance estimates for the 172 giant stars with 
I lo g gecheiie ~ lo g .Qrav e I < 0.5 are shorter than those 
of IZwitter et all (j2010) by 15% ± 24%, while the 73 
MS/SG stars have distances shorter by 36% ±21%. Note 
that our technique was optimized for metal-poor stars 
and uses parameter s derived from echelle spectra, while 
IZwitter et al.l (2010) optimized their method for all stars 
in the RAVE catalog, which have a high mean metallic- 
ity and typically younger ages and used parameter values 
from the RAVE pipeline analysis. 

5.6. Final Distances 

Our final estimate of the distance to each star, from the 
Sun, can be found in TableEl with the error being 20% of 
the distance estimate. The average distance to the RGB 
and RC/HB stars is ~ 2 kpc. All had distances less than 
~ 7 kpc, except one at ~ 16 kpc (see Figure fl"3|) . The 
MS/SG stars in our sample have an average distance of 
~ 220 pc, extending out to ~ 400 pc. The majority of 
our stars (primarily giants) extend to an average vertical 
height of \z\ ^ 1 kpc, as shown in Figure ITUl It is clear 
that our sample probes distances much further than the 
solar neighborhood (~ 100 pc) - not done by any previ- 
ous analyses of elemental abundances of metal-poor thick 
disk stars in the literature. 
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Figure 12. Plot of log g vs. log(T cff ) for the MS/SG stars in our 
metal-poor sample (after all parameter corrections) for different 
metallicity bins as shown in eac h p anel. The isochrones and bins 
are the same as those in Figure [TT] Stars with [Fc/H] < —0.8 are 
all old, except for two younger stars in the [Fe/H] = ( — 1,-0.8) 
bin. We therefore assume a mass of 1.0 Mq for the two younger 
stars, while the remaining stars with [Fe/H] < —0.8 are assumed 
to have a mass of 0.8 Mq. The ages of stars with [Fe/H] > —0.8 
are unclear, and cannot be easily separated. We therefore adopt a 
mass estimate of 0.9 Mq for these stars. 

5.7. Velocities and Orbits 

We computed three-dimensional space motions of our 
stars in cylindrical coordinates (given in Table[3]) by com- 
bining the distances and radial velocities derived from the 
analysis of our spectra with the proper motions given in 
the RAVE database. We further derived the full orbit of 
each star by assuming a Galactic potential and combin- 
ing that with the velocity and position information for 
each star. Stellar orbits were computed over 15 Gyr using 
an orbital integrator based on a three-component Galac- 
tic potenti al. The disk is modeled a s a Miyamoto-Nagai 
potential (iMivamoto fc Nagailli"975T ). while we used the 
iHernquistl (|1990f ) potential for the bulge. Finally, we as- 
sumed a logarithmic spherical potential for the halo. We 
took M disk = 8.0 x 10 10 M , M bulgc = 2.5 x 10 10 M , 
with characteristic scales a — 6.5, b — 0.26, c = 0.7, and 
al = 12.0, all in kiloparsecs, and w^ alo = 27000 (km s -1 ) 2 . 
These values ensure that the circular velocity equals 
220 km s _1 at 8 kpc from the Galactic center. 

The orbital parameters of each star are listed in Ta- 
ble [U Two sets of data are listed. The maximum 
vertical distance, |z ma £c|, and the closest and furthest 
distances, r pcr and r apo , reached by a star for all orbits 
integrated are listed in columns 3 through 5, while the 
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Figure 13. Distribution of total distance (upper panel) and ver- 
tical |z|-height (lower panel) for our candidate sample. Note that 
the star at a distance of ~ 16 kpc is not shown. 

same parameters for the final orbit of the star are listed 
in columns 6 through 8. The total number of orbits in- 
tegrated is given as iV or bit- The eccentricity of any given 
orbit of a star is defined as 

(14) 



The parameters for the last orbit were chosen to compute 
e in order to make direct comparisons of stellar eccentric- 
ities with model simulations (see The amplitude of 



variation for 



per 



and 



apo 



depends on the z-excursions 



during each orbit of a star, and thus, can be quite large. 
The amplitude of variation in the eccentricity of an orbit, 
however, is typically less than 0.05. We therefore assume 
the eccentricity of the last orbit is representative of the 
true orbital eccentricity of a star. 



6. POPULATION ASSIGNMENTS 

Each star was assigned to a Galactic population follow- 
ing the same Monte-Carlo method as in RIO. We drew 
10,000 random samples of each component of a star's 
space motion from a Gaussian error distribution centered 
on our estimate of the component velocity. The probabil- 
ities that each random set of velocities was drawn from 
the thin disk, thick disk, or halo were then computed us- 
ing the local characteristic Gaussian definitions for each 
Galactic population (given in Table 1 of R10 and repro- 
duced here in Table [SJ . A random set was assigned to 
a specific Galactic population if the probability was four 
times that of the other two probabilities. A star was then 
finally assigned to the Galactic population with the most 
random set assignments. 

It is important to note that some stars had probabili- 
ties that did not easily distinguish between the Galactic 
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Figure 14. Toomre diagram for our full sample with cyi < 
100 km s" 1 . The black circles, green diamonds, red squares, orange 
plus signs, and blue triangles correspond to thin disk, thin/thick, 
thick disk, thick/halo, and halo stars, respectively. The dashed 
curves indicate constant space motion, while an estimate of the lo- 
cal escape velocity, based on radial velocities within the RAVE 
database (Smith ct al. 200?5), is represented by the thick solid 
curve. Typical velocity errors are la error < 20 km s . 



populations (the ratio of the probabilities of two Galac- 
tic components was less than four). These stars were 
then assigned to additional intermediate thin/thick or 
thick/halo populations. This method, however, is sus- 
ceptible to possible misassignments. During the analy- 
sis, we therefore retained an additional probability statis- 
tic, which equalled the sum of the probabilities (normal- 
ized) obtained from all Monte-Carlo realizations for each 
Galactic population (hereafter, the PDF values). 

Similarly, random sets of a star's distance were drawn 
from a Gaussian error distribution centered on our esti- 
mate of the distance and a sigma equal to 20 percent of 
the distance. The probability that a random 'star' lies 
in each Galactic population was computed by comparing 
to the characteristic double-exponential distributions for 
the thin and thick disks and the two-axial p o wer-la w el- 
lipsoid for the halo, taken from Uuric et al.l ()2008[ ). As 
with the velocities, a star was assigned to the Galactic 
population with the highest number of occurrences from 
the re-sampled distance. The positional assignment was 
then used as a boundary condition such that if a star 
was assigned to the thick disk from its velocities and it 
was assigned to the halo from its position, the star would 
then be assigned to the halo. A star would remain as- 
signed to the thick disk, however, if it was assigned to 
the thin disk from its position. 

The stars were assigned as follows: 88 thick disk, 21 
thin disk, 51 intermediate thin/thick disk, 36 interme- 
diate thick/halo, and 123 halo. The last four columns 
of Table [3] gives the PDF values and Galactic population 
assignment for each star. A Toomre diagram is plotted in 
Figure [Ml to illustrate the relationship between our final 
population assignments and stellar velocities. Note that 
comparisons with the Toomre diagram from R10 show 
that the MS/SG stars mainly comprise the thin disk and 
thin/thick population. We further sum the stars' PDF 
values for each Galactic population within a given veloc- 
ity bin and plot the distribution of each velocity compo- 
nent in Figure [151 It is clear that the distributions reflect 
the underlying assumed populations. 



12 Ruchti et al. 



Table 3 

Kinematic Data 



Star 


D 


Vhelio 


Vn 


en 


v@ 




V z 




^thin 


-Pthick 


^halo 


POP 




(pc) 


(km s- 1 ) 


(km s- 1 ) 


(km s- 1 ) 


(km s 1 ) 


(km s- 1 ) 


(km s" 1 ) 


(km s- 1 ) 










C0023306-163143 


921 


-7.0 


53.2 


19.3 


-195.3 


85.9 


-77.9 


19.0 


0.00 


0.00 


1.00 


3 


C0315358-094743 


2434 


131.6 


86.3 


14.6 


50.4 


36.4 


-54.0 


15.2 


0.00 


0.13 


0.87 


3 


C0408404-462531 


15699 


52.3 


74.2 


344.8 


-88.9 


338.4 


241.7 


269.9 


0.00 


0.01 


0.99 


3 


C0549576-334007 


1157 


79.7 


-29.1 


16.4 


134.6 


11.7 


-24.7 


14.2 


0.02 


0.91 


0.07 


2 


C1141088-453528 


6152 


83.9 


-105.0 


54.2 


58.1 


51.9 


-88.1 


69.2 


0.00 


0.11 


0.89 


3 



Note. — Tabic [3] is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding 
its form and content. 



Table 4 

Orbital Parameters 









All Orbits 




Final Orbit 




Star 


Norbit 


''"per 




l^max | 


fper 


^*apo 


| ^max | 


e 






(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 




C0023306-163143 


100 


6.5 


9.2 


1.6 


6.6 


9.1 


1.5 


0.2 


C0315358-094743 


116 


1.4 


10.6 


2.3 


1.6 


10.5 


2.2 


0.7 


C0408404-462531 


61 


15.9 


23.8 


23.8 


16.9 


23.5 


23.4 


0.2 


C0549576-334007 


119 


3.9 


8.9 


0.6 


3.9 


8.9 


0.6 


0.4 


C1141088-453528 


127 


0.6 


10.6 


3.4 


0.7 


9.8 


2.4 


0.9 



Note. — Tabic ^ is published in its entirety in the electronic edition of the Astro- 
physical Journal. A portion is shown here for guidance regarding its form and content. 



Table 5 

Local Characteristic Velocity Distributions 



Population o~u o~q o~z Ref. 

(km s" 1 ) (km s" 1 ) (km s —1 ) (km s" 1 ) 



Thin Disk 39 20 20 -15 Soubiran ct al. (2003) 

Thick Disk 63 39 39 -51 Soubiran et"aT (2003) 

Halo 141 106 94 -220 Chiba fc B eers (2000) 



7. THE METAL-POOR THICK DISK 

7.1. Iron Abundance Distribution and Gradients 

Recall (from Q that we estimate the iron abundance 
of each star from Fen, since Fell is both the domi- 
nant species and is much less sensitive to non-LTE ef- 
fects than is Fel dTnevenin fc Idiart]|1999tlAsplund et al.l 
H999tlMashonkina et al.ll2011l ). FigurelT6lshows the PDF 
values for each Galactic population versus [Fe/H], while 
Figure [17] shows the distribution of iron abundance for 
each of the three Galactic populations. As was found in 
R10, the majority of the thick disk has [Fe/H] > —1.8, 
with a small tail to much lower metallicities, the lowest 
metallicity being [Fe/H] 2.7. The MS/SG stars as- 
signed to the thick disk predominantly have metallicities 
close to -1 dex, with a short tail down to [Fe/H] ~ —1.8. 
The distribution of [Fe/H], therefore did not show a sig- 
nificant change from that of R10. Due to sample selection 
effects, we do not completely sample the high-metallicity 
([Fe/H] > —1) parts of the distributions. Further, the 
offset between our final [Fe/H] values and the RAVE 
metallicities (see Figure fTU|) contributed to the shortage 
of thick disk stars with metallicities at the peak from lo- 
cal studies ([Fe/H] ~ —0.6). Stars that may have seemed 
to have metallicities near -0.6 dex, would end up with a 
lower value for the final metallicities. 



The most notable attribute of the thick disk metallic- 
ity distribution is that it appears to be double-peaked. 
The second, low-metallicity peak is comprised of many 
stars that have intermediate probabilities for both the 
thick disk and halo (PDFs < 0.6, assigned as thick/halo), 
while the peak above [Fe/H] ~ — 1 is comprised of stars 
with probabilities of being thick disk equal to 0.6 — 1.0 
(see Figure [T6|). It is therefore important to note that the 
low-metallicity peak is not just systematic noise in the 
population assignments, but is actually different. Could 
this be a metal-poor 'sub-component' of the thick disk, 
or could it be the high angular momentum tail of the 
halo? We will come back to this later (see §7.3p . 

Before we can determine the presence of iron abun- 
dance gradients in the thick disk, we must first determine 
possible selection effects in our sample and introduced 
through the analysis. For example, we introduced a tem- 
perature cut, such that only stars with temperatures in 
the range 4000 — 6500 K were selected for observation, to 
reduce contamination of stars that would most likely fail 
our abundance analysis (see . Figure fTHl plots [Fe/H] 
versus R and \z\ for our sample on top of the RAVE cat- 
alog, with [M/H] < —0.5, from which our sample stars 
were selected. We performed a linear fit to the [Fe/H] vs. 
[Fe/H] — [M/H] relation in Figure ITUl in order to put our 
sample stars and the RAVE catalog stars on the same 



Properties of the Thick Disk 



13 




-400 



-300 



-200 -100 
V n (km S" 1 ) 



100 



200 



300 



0.014 
0.012 
£0.010 
^0.008 
o 0.006 
£ 0.004 
0.002 
0.000 



-400 



-300 



-200 -100 
V @ (km S~ 1 ) 



100 



200 



300 




-200 -100 
V z (km s~ 1 ) 



300 



Figure 15. Distribution of Vn (top panel), Ve (middle panel), 
and Vz (bottom panel) velocities for our stars given the PDF value 
that they belong to the thin disk (solid black), thick disk (dashed 
red), and halo (dot-dashed blue). Note that each distribution was 
created by summing the PDF values within a given velocity bin. 
The distributions were then normalized such that the total area 
equals unity. 

scale. Distances to the RAVE catalog stars were com- 
puted using an analogous method to that described in 
Our sample does not show the same distance coverage as 
the RAVE catalog stars for [Fe/H] > -1.2. This implies 
that the lack of thick-disk stars with [Fe/H] > —1.2 at 
large radial and vertical distances is possibly an artifact 
of our cuts in temperature. 

We investigated the possible origin of these selection 
effects by comparing to an old (12 Gyr) isochrone. In 
Figure [T51 we plot the derived data for our sample to- 
gether with 12 Gyr isochrones of differing metallicities. 
The metallicities of the isochrones were simply converted 
to [Fe/H] by [Fe/H] = \og(Z/Z Q ) as a first-order ap- 
proximation. We applied an apparent magnitude limit 
of I = 10, the peak /-magnitude of our sample (see Fig- 
ure [3]) , to the isochrone data to determine the maximum 
distance that can be observed for each metallicity. It is 
clear that the stars of our sample do not reach the max- 
imum limit (as shown in the top panel of Figure \T§\) . 
However, when we apply the same temperature cut that 
we applied during our candidate selection (4000—6500 K) 
to the isochrone data, then we see a similar trend emerge 
in the isochrone behavior as in the sample (lower panel, 
Figure [TH]). This is a clear illustration of the selection 
effects in our final sample. For [Fe/H] < —1.2, however, 
the isochrone data are unaffected by the temperature cut, 
suggesting that below this value the data for our sample 
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Figure 16. PDF values for each Galactic population versus 
[Fe/H]. 

may be used without further corrections. 

Only the 49 stars assigned to the thick disk with 
[Fe/H] < —1.2 were used to assess the amplitude of 
metallicity gradients in our metal-poor thick-disk sam- 
ple, as shown in Figure [20] This figure shows the robust 
least squares fits to the data along with 95% (2 sigma) 
confidence intervals for each fit. Both fits have slopes 
which are formally non-zero, but not significant, corre- 
sponding to a —0.09 ± 0.05 dex kpc" 1 gradient in the 
vertical direction and a +0.01 ± 0.04 dex kpc -1 gradient 
in the radial direction. Iron abundance versus |z ma x| an d 
fapo (values are maxima for all integrated orbits) also 
exhibits similar gradients (see Figure [2~Tj) . 

As stated previously, the population assignments are 
susceptible to possible misassignments. We therefore 
checked our results using our PDF values. The metal- 
licity with the maximum sum of thick-disk PDFs inside 
a specific velocity bin was chosen as the metallicity in 
that given bin. We then fit a line across the maxima to 
determine possible gradients. The results showed very 
similar results are very similar to those found using only 
stars assigned to the thick disk. 

It is important to note that at low iron abundance and 
large distances we are hindered by small-number statis- 
tics for the thick disk. We therefore set up a bootstrap 
analysis, in which we created 10,000 re-samples consist- 
ing of 25 stars randomly selected from the 49 stars as- 
signed to thick disk with [Fe/H] < —1.2. We performed 
a least squares fit to each re-sample, and then took the 
average and standard deviation of all resamples to de- 
termine the range to which the fit is affected by possible 
outliers. 
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-2.0 
[Fe/H] 

Figure 17. Distribution of [Fe/H] for thin disk, thick disk, and 
halo. Note that each distribution was created by summing the PDF 
values within a given metallicity bin. The distributions were then 
normalized such that the total area equals unity. Recall that our 
sample was selected to be metal-poor, so that the full metallicity 
distributions of each component are not uniformly represented. 

Figure [22] shows the results from this test. The con- 
fidence intervals shown in this figure represent the de- 
gree to which the fit changes for different random sam- 
ples. The mean slopes are now only slightly steeper 
than for our original fits, giving a vertical gradient of 
-0.13 ± 0.07 dex kpc -1 and a +0.06 ± 0.06 dex kpc -1 
radial gradient. Similarly to our original fits, a slope 
equal to zero is also still consistent within errors. 

7.2. Alpha-to-Iron Ratios 

Figure [23] displays several [a/Fe] ratios versus [Fe/H] 
for our entire sample, which can be directly compared to 
the same plots in R10 for only the giant stars. Most of 
the stars with [Fe/H] > -0.7 consist of MS/SG stars that 
were assigned to either the thin disk or thin/thick pop- 
ulation, and have an a-enhancement lower than for the 
stars at lower metallicities. The more metal-poor MS/SG 
stars, however, typically have similar a-enhancement to 
the giants of the same metallicity, an indication that the 
results from R10 were not affected by the addition of the 
MS/SG stars. 

The metal-poor thick-disk stars have [a/Fe] signif- 
icantly above solar, ~ 0.36 dex for Mg and Si and 
- 0.27 dex for Ca and Tin. The mean [Til/Fe] value 
is about a tenth of a dex lower than that for Ti II. This 
offset is most lik ely due to non-LT E effects present within 
our analysis f see iBergem ann 201l|). The ratios also show 
low scatter, ;$ 0.09 dex for all a-elements, which is less 
than the 0.1 dex experimental error in [a/Fe]. The [a/Fe] 
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Figure 18. [Fe/H] vs. Ft and |z| for our sample of stars (symbols 
and colors are the same as Figure [141) . The small data points are 
those stars from our original RAVE catalog, from which our sample 
was taken. The RAVE catalog stars populate large |zj, even at 
high [Fe/H], which suggests that the lack of stars in our metal- 
poor sample with [Fe/H] > —1 at large |z| is due to the selection 
function of the sample, combined with the subsequent temperature 
cuts prior to the abundance analysis. 



ratios also blend smoothly into the halo stars with the 
difference in mean [a/Fe] ranging between 0.00—0.03 dex, 
well within our experimental error. There are a few 
thick-disk stars with [Fe/H] > —0.7 that may have lower 
enhancement than the metal-poor thick-disk stars, but 
they do not represent a large enough sample to make any 
clear conclusion. Also note that, as was shown in R10, 
there is a thick disk star at [Fe/H] ~ — 1 with very large 
[Si/Fe] enhancement and a halo star at [Fe/H] ~ —1.45 
that shows consistently low a-enhancement. These stars, 
however, show no peculiarities in their kinematics, and 
will be the subject of future papers. 

Figures [Ml and [231 show the [a/Fe] ratios versus R and 
\z\, respectively, for only stars with [Fe/H] < —1.2. In 
both plots, the stars assigned to the thick disk show lit- 
tle dependence on position, for all a-elements. In the 
vertical direction, [Mg/Fe] and [Si/Fe] slightly increase 
towards larger \z\ at +0.03±0.02 dex kpc" 1 and +0.02± 
0.01 dex kpc -1 , respectively, while [Ca/Fe] decreases at a 
rate of -0.01±0.01 dex kpc -1 and [Til/Fe] and [Tin/Fe] 
decrease at a rate of —0.02 ± 0.02 dex kpc -1 . The 
ratio of all five elements to iron increase at less than 
0.03 ±0.01 dex kpc -1 radially outward. As in the case of 
iron abundance, the [a/Fe] ratios versus r^po and |z ma x| 
show similar gradients (see Figures [231 and [27]) . only shift- 
ing by H 0.01 dex kpc -1 vertically. In addition, the 
radial gradient reduced to ~ 0.01 dex kpc -1 for all a- 
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Figure 19. [Fc/H] vs. Distance for our sample of stars. The red 
'plus' symbols represent a 12 Gyr isochrone of varying metallicity, 
for which distances are computed assuming an apparent magnitude 
of J = 10. The black points represent our sample of candidate 
metal-poor thick disk stars. Note that there is one star not shown 
at a distance of 15 kpc with [Fe/H] ~ —2.2. The top panel shows 
the isochrone points without any temperature cut, while the same 
temperature cut as that of our candidate sample (4000 — 6500 K) is 
applied to the isochrone data in the bottom panel. The brightest 
metal-rich giants are cooler than their metal-poor counterparts, 
and so the low-temperature limit cuts them out and makes the 
difference in isochrones seen in the two panels. For [Fe/H] < — 1.2 
(blue line) the sample is unaffected by our temperature cut. 

elements. 

7.3. High Angular Momentum Halo 

In Figure it is evident that some stars assigned to 
the thick disk have overlapping Vq velocities with the 
high-velocity tail of the halo. Further, we found that 
there is a second, low-metallicity peak in the metallicity 
distribution of the thick disk (Figure |T7| . 

We investigated the possibility that the low metallic- 
ity peak might indicate contamination from the high- 
angular momentum halo by plotting [Fe/H] and each 
of the [a/Fe] ratios versus Vg>, which are shown in Fig- 
ure [28j Within the thick disk, there is no trend of [a/Fe] 
or [Fe/H] with Ve, except perhaps at the regime when the 
thick disk overlaps with the thin disk. This implies that 
there is no difference between those stars that might kine- 
matically be a part of the tail of the halo and those that 
have azimuthal velocities that are too high to be halo. 
Thus, Figure [28] shows there is no difference between the 
halo and thick disk, as was found when comparing the 
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Figure 20. [Fe/H] vs. R and \z\ for stars only assigned to the 
thick disk (with [Fe/H] < —1.2) that are not biased by our se- 
lection effects, represented as the black points. Typical errors are 
~ 0.1 kpc in R and |z| and 0.1 dex in the [a/Fe] ratios. The black 
line is the robust least squares fit to the data, while the red-dashed 
curves represent the 95-percentile confidence intervals for the fit. 
Note that a slope equal to zero (no gradient) is consistent with our 
data for both R and \z\. 



-1.0 
-1.5 


- 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

. ', S% o ° o o 


1 ' ' ' 


i I i i i i I i i i. 


£-2.0 






-° — 




-2.5 










-3.0 


- . i . , 







r 




0.5 


1.0 1.5 2.0 

IZmaxI (kpc) 


2.5 


3.0 3.5 


-1.0 


1111,1 


i i i | i i i i | i i i i | i i i 


> 1 > 


i I i i i i I i i i i 














; — as 


Z 3 - -** -it - - 




6 ■ 

a '. 


£-2.0 










-2.5 










-3.0 


- 







~ 




6 


7 8 9 


10 


11 12 



r ; , po (kpc) 

Figure 21. [Fe/H] vs. r apo and |z m ax| (columns 4 and 5 of Ta- 
ble [4j for stars assigned to the thick disk with [Fe/H] < —1.2. 
Typical errors are ~ 0.1 kpc in r apo and |z max | and 0.1 dex in 
the [a/Fe] ratios. The black line is the robust least squares fit to 
the data, while the red-dashed curves represent the 95-perccntile 
confidence intervals for the fit. 

[a/Fe] ratios versus [Fe/H]. 

In Figure [29] we plot the velocity distributions for only 
stars with [Fe/H] < —1.2. The thick disk distributions 
still show a clear distinction from the halo distributions. 
It is therefore reasonable to assume these stars could still 
be thick disk stars. 

7.4. A Lagged Thick Disk Component? 

Some have proposed (e.g., iCarollo et~ al. 2010) that 
the kinematics of the metal-poor thick disk may dif- 
fer from the canonical thick disk. In this case the 
metal-poor component would have higher velocity dis- 
persions and a slower rotational velocity. To investi- 
gate the effects of such a component on our results, we 
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Figure 22. [Fe/H] vs. R and |z| for stars assigned to the thick 
disk (with [Fe/H] < —1.2). The black lines are now the mean 
of least squares fits to 10,000 random samples, consisting of only 
half of the stars assigned to the thick disk. The red-dashed curves 
represent the 95% confidence intervals derived from the scatter in 
the distribution of linear fits to each random sample. 
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Figure 23. Computed [a/Fe] ratios vs. [Fe/H] for our sam- 
ple. Note that element ratios are computed using the iron abun- 
dance of the same ionization st ate as the a-el e ment (e.g., [Si/Fe] 
= [Sil/Fel]) as is suggested by Kraf t & Ivans (2003). Color and 
symbols are the same as in Fig. 1 1 41 The cause of the offset be- 
tween [Til/Fe] and [Tio/Fe] is most likely non-LTE effects present 
within the abundance analysis, but the thick disk and halo still 
show similar enhancement in each. 

adopted a (Ve) = —100 km s _1 and velocity dispersions 
of (crn, ^e, Qz) = (63, 70, 60) km s -1 for the lagged com- 
ponent fsee lGilmore et aill2002D in addition to the thin 
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Figure 24. Computed [a/Fe] ratios vs. R for our sample of stars 
with [Fe/H] < —1.2. The left plots s how all populations with 
symbols and colors the same as Figure 1231 The right plots show 
the least squares fits to only stars assigned to the thick disk. The 
black line is the fit, while the red-dashed curves represent the 95% 
confidence intervals. 



disk, canonical thick disk, and halo components. We 
then computed new population assignments for the stars. 
This increased the number of stars assigned to the thick 
disk population by 22 stars; 16 of which were formerly 
assigned to the thick/halo population, 4 formerly halo, 
and 2 formerly thin/thick. 

The overall shape of the iron distribution did not 
change for the combined thick disk components. Further, 
the mean [a/Fe] ratios showed no difference from that of 
the canonical thick disk alone, while the gradients in both 
[Fe/H] and [a/Fe] changed by less than 0.01 dex kpc" 1 . 
Our results therefore do not show any significant change 
by including a lagged component for the thick disk. 

8. IMF VARIATION 

In R10, we established that our metal-poor thick-disk 
giants stars formed during a period of rapid star forma- 
tion, primarily pre-enriched by core-collapse supernova 
(e.g. SNe II). In this section we quantify the level of 
IMF variation (specifically the slope, x, of the IMF) of 
the massive stars that ended as a core-collapse SNe using 
model yields and comparing with the scatter in our data. 

For this test, we adop ted the mass-dependent Mg and 
Fe yields of SNe II from lKobavashi et all (|2006l ). we can 
then compute the massive-star IMF-averaged yield for 
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Figure 25. Computed [ct/Fe] ratios vs. z for our sample of s tars 
with [Fc/H] < —1.2. Symbols and colors are the same as Figure l24l 
Recall that the left-hand side shows all populations, while the right- 
hand side shows only the stars assigned to the thick disk. 

each element using different IMF slopes, a;, by: 



Yimf = 



Mo 



Y(M) M~ (1+x) dM 



(15) 



where Y{M) is the mass-dependent yield, M\ — 13 M , 
and M up = 40 M Q . 

Figure [301 shows a plot of [Mg/Fe]iMF versus the IMF 
slope, x. We used this plot to determine the spread in 
IMF slope values from the scatter in the [Mg/Fe] ratios 
of the stars assigned to the thick disk. 

From Figure [SHI a small difference in IMF slope implies 
a large difference in [a/Fe] (in this case [Mg/Fe] imf), bet- 
ter than scatter. Previously, we showed that the thick 
disk has no scatter outside of random errors in all [a/Fe] 
ratios. This leaves no room for any variation in the 
IMF. Further, this implies that the ISM was well-mixed 
prior to star formation. Similarly, the difference between 
the mean [Mg/Fe] values of the halo and thick disk is 
0.03 dex, well within our 0.1 dex errors. This illustrates 
that the halo and metal-poor thick disk came from a very 
similar massive star IMF. 

9. ORBITAL ECCENTRICITY 

I Sales et "all ()2009L hereafter S09) investigated the util- 
ity of orbital eccentricity (e) distributions as a tool for 
distinguishing several scenarios for the formation of the 
thick disk. They compared the predictions from several 
simulations of the formation of thick disks and found 
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Figure 26. Computed [a/Fe] ratios VS. 7"apo (column 4 of Tablc[4j 
for our sample of star s wi th [Fe/H] < —1.2. Symbols and colors 
are the same as Figure [24] Recall that the left-hand side shows all 
populations, while the right-hand side shows only the stars assigned 
to the thick disk. 

that the e-distributions also provide a robust diagnos- 
tic to distinguish between stellar populations that form 
in the simulated galaxy (in situ) versus those that are 
formed in a satellite galaxy and accreted into the sim- 
ulated galaxy for all scenarios that involve both types 
of populations. The accreted population dominates high 
eccentricities, while the in situ population dominates the 
lower eccentricity bins. It is important to note that only 
one realization, with a specific set of initial conditions, 
for each formation scenario was used to compute the ec- 
centricity distributions. This is, however, sufficient since 
there is no reason that the given simulations are not rep- 
resentative of that scenario. 

It is important to note that we cannot use our pop- 
ulation assignments when constructing the distribution 
of orbital eccentricities for the thick disk. The shape 
and extent of the tail of the thick-disk distribution is 
biased by our definition of 'thick disk' during our pop- 
ulation analysis. We simulated our analysis by creating 
1000 model stars in which their distance and Galactic 
coordinates were randomly chosen assuming a uniform 
distribution, and their 3-dimensional space motion was 
randomly selected using the 'thick disk' Gaussian defini- 
tion in Table [5j Each model star was then run through 
our orbital program to determine the orbital eccentricity. 
Only 12% of these model thick disk stars had e > 0.6, 
which indicates that our definition of thick disk does not 
allow for many stars on highly eccentric orbits. Further, 
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Figure 27. Computed [a/Fe] ratios vs. |z max | (column 5 of Ta- 
ble [4} for our sample of stars wit h [Fe/H] < —1.2. Symbols and 
colors are the same as Figure 1241 Recall that the left-hand side 
shows all populations, while the right-hand side shows only the 
stars assigned to the thick disk. 



this could also affect the shape of the right side of the 
eccentricity distributions. 

We therefore followed the procedure of S09 to inves- 
tigate the orbital eccentricities of our metal-poor thick 
disk sample stars, calculated in iJ5.7l for the last orbit of 
a star. In this case, we computed the e-distribution of 
stars with V& > 50 km s _1 and 1 < \z\/ z& < 3, where 
z ( \ = 0.9 kpc. This distribution is shown in Figure [3T1 
To estimate the effects of distance errors on the distribu- 
tion, we also plot the resultant eccentricity distribution 
when distances are increased and decreased by 20%. 

The majority of the stars in this sample appear to 
exhibit low orbital eccentricities, with the distribution 
peaking around s ~ 0.3 — 0.4. Our distribution does 
not show a strong resemblance to any of the distribu- 
tions found in S09. The tail to high eccentricity some- 
what resembles th at of the direct accretion scenario of 
lAbadi et al.l ((2003). The peak at lower eccentricities is 
significantly lower than predicted by the direct accre- 
tion scenario, and is more consistent with the predic- 
tions of models wherein the stars of the thick disk were 
primarily formed in situ. It is possible that we are under- 
sampling the low-eccentricity bins since we do not have 
many stars at the peak of the thick-disk metallicity dis- 
tribution, [Fe/H] ~ — 0.6 . It is likely that including more 
metal-rich thick-disk stars could significantly change the 
distribution of eccentricities. 
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Figure 28. [Fe/H] and [a/Fe] ratios vs. V© for ou r sample of 
stars. Symbols and colors are the same as Figure l23l Note that 
the error in V© is typically lower than the error bar shown. 

10. DISCUSSION 

In this work, we analyzed our full sample of metal-poor 
thick disk candidate stars. The MS/SG stars, added 
to our giant star sample from R10, were primarily as- 
signed to the thin disk and thin/thick populations, and 
increased the number of stars with [Fe/H] > —0.8. We 
found that our findings from R10 were unchanged, in 
that the [a/Fe] ratios for the metal-poor thick-disk stars 
are enhanced and show low scatter (< 0.09 dex, within 
the error of 0.1 dex), indicating that star formation took 
place on a short timescale in which the metal-poor thick 
disk was pre-enriched by core-collapse SNe from an in- 
variant massive star IMF. Further, the metal-poor thick 
disk and halo were most likely pre-enriched by the same 
massive star IMF, showing a difference in a-enhancement 
of < 0.03 dex. The low amplitude of scatter in the ele- 
ment ratios indicates that the ISM from which the stars 
formed was well-mixed. 

As discussed in R10, the enhancement and low scatter 
of [a/Fe] in the thick disk are evidence that the formation 
of the thick disk had little influence by the late, direct ac- 
cretion of stars from dwarf galaxies. The a-enhancement 
in the metal-poor thick disk contrasts with the expecta- 
tions from models that h ave direct accretio n up until 
about ~ 6 Gyr ago (e.g.. lAbadi et all 12003:. assuming 
that the accreted dwarfs formed stars and self-enriched 
similarly to the surviving dwarfs. In this case, the stars 
accreted would then look chemically different from our 
metal-poor thick-disk stars (see R10). Further, the dis- 
tribution of orbital eccentricities for our stars does not 
resemble that of the direct accretion scenario, instead 
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Figure 29. Distribution of Vn (top panel), Ve (middle panel), 
and Vz (bottom panel) velocities for stars with [Fe/H] < —1.2, 
given the PDF value that they belong to the thin disk (solid black), 
thick disk (dashed red), and halo (dot-dashed blue). Note that 
each distribution was created by summing the PDF values within 
a given velocity bin. The distributions were then normalized such 
that the total area equals unity. 
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Figure 30. IMF-aver aged [Mg/Fe]iMP vs. the IMF slope, x, com- 
puted using equation 1151 . Note that as x decreases, the ratio of 
Mg-to-Fe increases, and that the relationship is nearly linear. The 
Salpeter IMF occurs for x = 2.35. 

resembling a population that was formed i n situ. 

Usin g the full RAVE catalog of stars, iWilson et al.l 
(|2011D also computed orbital eccentricities for a more 
uniform sample of thick disk stars, finding results consis - 
tent with ours. On the other hand. lDlerickx et al.l (|2010f ) 
looked at the orbital eccentricities of a sample of thick- 
disk candidate stars selected from the Sloan Digital Sky 
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Figure 31. The distribution of orbital eccentricities for those stars 
in our sample with Vg> > 50 km s _1 and 1 < |z|/«d ^ 3, where 
Zd = 0.9 kpc. The distribution when the distances to the stars are 
decreased or increased by 20% are given by red, dashed and blue, 
dot-dashed lines, respectively. 

Survey, and found that their distribution is inconsistent 
with the thick disk forming from radial migration only 
or hea ting due to mergers only. iCasetti-Dinescu et al.1 
(|2011|) . however, show that their distribution, using data 
from RAVE and the Southern Proper Motion Program, 
is most consistent with the heating and merger scenar- 
ios. Overall, these results still support our finding that 
the stars in the thick disk primarily formed in situ. It 
is also important to note that our distribution of eccen- 
tricities is the first to have the metal-poor thick disk 
well-represented, but selection biases may be affecting 
the shape of the distribution. 

From previous investigations of the local neighborhood 
and our more extended sample, it is clear that the thick 
disk is old, and had to form stars during a short, rapid 
burst to produce the a-enhancement we detect. Direct 
accretion of stars from dwarf galaxies that formed a long 
time ago (about 10 Gyr ago) during a period of rapid 
star formation is then still viable. These types of dwarf 
galaxies, however, are extremely rare. All known dwarf 
spheroidal galaxies typically have extended star forma- 
tion and lower a-enhancement at the sam e [Fe/H] range 
of the metal-poor thick disk (see review bv lTolstov et al.l 
2009). We can therefore conclude with confidence that 
the late accretion of stars from satellite galaxies did not 
play a role in the formation of the thick disk. 

Models that include a significant component of the 
thick disk formed in situ, however, need to be further 
assessed. These models all imply or directly predict a 
thick disk with high a-enhancement. Further, given the 
errors in the distances to our stars, and the effects of 
our definition of 'thick disk', the orbital eccentricity dis- 
tribution of our sample of metal-poor thick disk stars 
does not show any direct inconsistencies with the distri- 
butions in S09 for these scenarios. It is interesting to 
note, however, that the S09 e-distribution for the heat- 
ing scenario shows a significant fraction of accreted stars 
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that comprise the thick disk. Unless all of these accreted 
stars were assimilated into the thick disk at early times 
(as stated above), then this fraction of stars would have 
lower a-enhancement than we see in our sample. 

The models also show differences in their predictions 
for radial or vertical abundance variation. Observational 
evidence for abundance gradients in the thick disk are 
therefore very important. The majority of previous stud- 
ies with high resolution analyses have shown no evidence 
for a vertical metallic i ty gradient in the thick d isk (cf., 
iMishenina et al J 120041: iSoubiran fc Girardll2005[ ). These 
studies, however, primarily include nearby stars with 
[Fe/H] > —1.5, while our gradient is computed for thick- 
disk stars with [Fe/H] < —1.2. More recent studies have 
shown evidence for the possibility of a vertical metallicity 
gradient. Looking at stars at high Galactic latitudes and 
z < 4 kpc. lSiegel et all (|2009f ) found that iron abundance 
decreases with vertical height by —0.15 dex kpc" 1 . This 
data, how ever, may have str ong contamination from the 
thin disk. Ilvezic et all (|2008l ) found a 0.1-0.2 dex kpc -1 
metallicity gradient for stars selected from SDSS that lie 
a vertical heights, z = 1 — 2 kpc. At these heights, the 
likelihood of contamination from thin disk is small, and 
so one could deduce that they were looking at mostly 
thick disk stars. 

The stars in our sample assigned to the thick disk lie 
within about ±2 kpc radially about the Sun, and have 
vertical distances of \z\ < 3 kpc. Due to biases intro- 
duced by our analysis technique, we only investigated 
abundance gradients for those stars with [Fe/H] < —1.2. 
We found very small amplitude radial and vertical gradi- 
ents, < 0.03 ± 0.02 dex kpc -1 , in the a-enhancement of 
the metal-poor thick disk. This further verifies that the 
ISM was well mixed during star formation. Further, we 
found a small +0.01 ± 0.04 dex kpc" 1 radial iron abun- 
dance gradient, however, it is possible that [Fe/H] signif- 
icantly changes by ~ —0.09 ±0.05 dex kpc" 1 with verti- 
cal height above the Galactic plane. A fit resulting in no 
gradient is still possible, however, within 2-sigma confi- 
dence limits. Note that it is possible that including the 
more metal-rich thick disk would significantly change the 
amplitude of the slope of the gradients. Our results, how- 
ever, resemble those found by groups studying the more 
metal-rich thick disk (as given above). 

A vertical metallicity gradient is expected from the 
dissipational collapse model. Additionally, a thick disk 
with uniformly enhanced a-abundances is also likely if 
the collapse took se veral millions of years, as predicted 
(|Burkert et al.|[T992l) . Rapid star formation early on dur- 
ing the heating scenario will produce enhanced [a/Fe] 
ratios, while small amplitude radial abundance gradients 
are also possible. The merger model also predicts a high 
star formation rate from gas dissipation, and hence en- 
hanced [a/Fe] ratios, but the final thick disk is expected 
to have uniform abundance ratios. Although our data 
does not directly disagree with a uniform disk, a signif- 
icant iron abundance gradient would challenge this sce- 
nario. It is also important to note that there is direct 
accretion of stars during the merger and heating scenar- 
ios, which can occur until late times. Late accretion of 
stars must therefore be extremely minimal for the abun- 
dances of the final thick disks in these scenarios to match 
the low scatter in [a/Fe] for our thick disk stars. 



An a-enhanced thick disk is also implied from the 
clumpy disk model in which there is rapid star forma- 
tion from internal processes with rapid mixing by strong 
turbulence. Additionally, the thick disk in the radial mi- 
gration model has high [a/Fe] ratios, since the stars in 
the inner disk und ergo rapid star formation a nd then 
move outward, (see iSchonrich fc Binnevl[2009bf ). No ra- 
dial abundance gradients are predicted due to blurring 
across radii in the radial migration model and the tur- 
bulent mixing in the clumpy disk model. Note that the 
radial migration model does not include effects from the 
bar of the Galaxy. It is possible that the combination 
of the bar and spiral arms could result in a very effi- 
cient mixing mechanism, which could cause variation in 
the metallicity of the thick disk as a function of radius 
(jMinchev fc Famaevl2010HMinchev et al.ll2011| ). It is un- 
clear, however, if a vertical iron abundance gradient is 
possible in either the radial migration model or clumpy 
disk model, indicating that more modeling must be done 
to ascertain this chemical signature. 

11. CONCLUSIONS 

The metal-poor thick disk of the Milky Way Galaxy 
is enhanced in the a-elements and reaches to metallici- 
ties down to -2 dex. We find that the stars in the thick 
disk most likely formed within the potential well of the 
Milky Way Galaxy. Direct accretion of stars could have 
occurred at very early times (~ 1 Gyr after the start 
of star formation) in the formation of the thick disk, but 
the later contribution of accreted stars into the thick disk 
was very minimal. The abundance trends of the metal- 
poor thick disk tends to favor models which result in a 
thick disk with a significant vertical metallicity gradient, 
however, a uniformly enhanced thick disk is still possi- 
ble. Additionally, the radial migration and internal star 
formation at high redshift may have also contributed to 
the formation of the thick disk, but more work needs 
to be done to quantify the abundance trends for these 
scenarios. 
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